Transport of fractional Hall quasiparticles through an antidot 
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Current statistics of an antidot in the fractional quantum Hall regime is studied for Laughlin's 
series. The chiral Luttinger liquid picture of edge states with a renormalized interaction exponent 
g is adopted. Several peculiar features are found in the sequential tunneling regime. On one side, 
current displays negative differential conductance and double-peak structures when g < 1. On 
the other side, universal sub-poissonian transport regimes are identified through an analysis of 
higher current moments. A comparison between Fano factor and skewness is proposed in order 
to clearly distinguish the charge of the carriers, regardless of possible non-universal interaction 
renormalizations. Super-poissonian statistics is obtained in the shot limit for g < 1, and plasmonic 
effects due to the finite-size antidot are tracked. 

PACS numbers: 73.23.-b,72.70.+m,73.43.Jn 



I. INTRODUCTION 

The peculiar properties of quasiparticles (qp) in the 
fractional quantum Hall effect (FQHE) have received 
great attention especially for the states at filling factor 
v = 1/p, p odd integer, in which gapped bulk excitations 
were predicted to exist and to possess fractional charge 1 
e* = ve (e < electron charge) and statistics^ 

A boundary restriction of this theory was subsequently 
put forward in terms of edge states by Wen£ This theory 
recovered the fractional numbers of quasiparticles in the 
framework of chiral Luttinger Liquids (%LL), and indi- 
cated tunneling as an accessible tool to probe them4 Ac- 
cordingly, quasiparticles with charge e/p were measured 
in shot noise experiments with point contact geometries 
and edge-edge backscattering^ 

A key prediction of %LL theories is that the interac- 
tion parameter should be universal and equal to v. As a 
consequence, the quasiparticle (electron) local tunneling 
density of states obeys a power law in energy D oc E v ~ x 
(D oc E 1 / u ~ 1 ). Several geometries have been set up in ex- 
periments to test this nonlincarity through measurements 
of tunneling current / versus bias voltage V. For in- 
stance, in the case of electron tunneling between a metal 
and an edge at filling v, one should have / oc V a in the 
limit eV 3> ksT with a = \jvX Experiments^ at fill- 
ing factor 1/3 indeed proved a power-law behaviour but 
with tt/:3. Deviations were observed also with quasipar- 
ticle tunneling in an almost open point contact geometry 
at v = l/3jl here, the predicted backscattering linear 
conductance is Gb oc T 2v ~ 2 while the measured quan- 
tity obeys a power law with positive exponent. Anal- 
ogous discrepancies were observed in a similar geome- 
try, with the quasiparticle tunneling differential conduc- 
tance developing a minimum around zero bias instead 
of a maximum for decreasing temperature.— Moreover, 
several numerical calculations and simulations also dis- 
agree with the conventional chiral Luttinger theories. We 
mention for instance finite-size exact-diagonalization cal- 
culations with short-range 9 - or Coulomb electron-electron 



interaction s 10 : 11 ' 12 and Monte Carlo simulations with in- 
teracting Composite Fermions on a ring. 13 

The disagreements of xLL predictions with ob- 
served exponents are still not completely understood, 
although several theoretical mechanisms have been 
put forward to reproduce a renormalized Luttinger 
parameter, including coupling to phonons or dis- 
sipative environments ) 1 W 6 ' 17 effects of interaction 
range ) 18 ' 19 edge reconstruction with smooth confining 
potentials^ 2 ^ 2 ^ 3 - 

Our purpose is to discuss fractional Hall edges in an 
enriched X-LL theory where the possibility of a renormal- 
ized interaction parameter g ^ v is assumed, analysing 
different transport regimes and clearly distinguishing sig- 
natures of charge ve qp from effects due to the quasipar- 
ticle propagators governed by g. To do so, we choose a 
geometry where two fractional quantum Hall edges are 
connected via weak tunneling through an edge state en- 
circling an antidot as in Fig. [T} 24 : 25 Such a setup has 
proven to be extremely versatile and controllable. It has 
been for instance employed in a series of experiments 
where fractional charge and statistics of quasiparticles 
were addressed^ The same geometry has been also used 
to detect blockade effects and Kondo physics with spinful 
edges in the integer regime.— 

As mentioned before, a series of experimental obser- 
vations of fractional charge has been based on noise 
measurements^ the possibility to extract a significant 
e*/e from such measurements is given by the fact that the 
geometry of the setup and the tunneling regime ensure 
a poissonian process. Smooth evolution in the charge 
carrier at v = 1/3 from e* = ve to e was observed in 
point contacts changing the backscattering amplitude via 
gates from to 1j2£ This appears consistent with evolu- 
tion from fractional qp to electron tunneling. 29 Anyway, 
this conclusion can only be obtained under the additional 
hypothesis of independent particle tunneling. 30 Only the 
observation of higher moments, as e.g. the normalized 
skewness, could cross-check this hypothesis.— Therefore, 
to fully explain the charge measurements one has to ob- 
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FIG. 1: Geometry of the system. In white, the Hall fluid; 
in black, the depleted areas defining the two quasiparticle 
tunneling points between the left (L) or right (R) edge and 
the antidot (^4). In this setting, the magnetic field points out 
of the plane. 



serve higher moments beyond the Fano factor. Recent 
advances in measurement techniques could open this in- 
triguing possibility, especially in view of the unexpected 
results recently reported.—^ 

In this paper we propose to compare Fano factor and 
skewness in transport regimes where the statistics of the 
charge transfer is not poissonian. We focus our attention 
on transport through an antidot in weak-backscattering 
and sequential tunneling limit, at fractional filling factor 
v = 1/p (p odd integer). Our task is twofold: on one 
side, we analyse the tunneling current. It presents re- 
markable features driven by g, like e.g. negative differen- 
tial conductance and double-peak structures. Power-law 
behaviours exist and can be used to determine the renor- 
malized interaction parameter. On the other side, we de- 
rive a method to assess fractional charge independently of 
possibly renormalized g =/= v. We analyse noise and skew- 
ness in processes with different transport statistics both 
in the shot and in the thermal limit. We find universal 
points that unambiguously define the fractional charge. 
In addition, we describe transport regions where the Fano 
factor is sensitive to the power laws of the quasiparticle 
propagators and presents super-poissonian correlations. 

The paper is organized as follows: in Sec. II the model 
is introduced and bosonization procedures are briefly re- 
viewed. Sections III and IV contain the description of 
methods and results. Finally, in Sec. V we discuss 
our findings and comment existing experimental appli- 
cations. 



II. MODEL 

In our model, edge states form at the boundaries of 
the sample and around the antidot (Fig. [1]); mesoscopic 
effects are associated to the finite size of the antidot 
through an Aharonov-Bohm (AB) coupling, and tunnel- 
ing barriers couple the circular antidot with both edges. 



The complete Hamiltonian reads 

H = Hl + Hl + Hl + H^+Hi + Hl, (1) 
and the individual pieces are now described in detail. 

A. Bosonization of free Hamiltonians 

In Eq. ([1]), Hi are standard Wen's hydrodynamical 
Hamiltonians for the left, right and antidot edge (I = 
L,R,A): in terms of the electron excess density pi(x) 
(H= 1), one has^s^ 

wit r^-i/ 2 

Hf = - dxpf(x), (2) 

v J-C/2 

where v is the edge magnetoplasmon velocity and Ci is 
the edge length. The theory is bosonized with the pre- 
scription pi(x) = d x cj)i(x)/2iT where <fii{x) are scalar fields 
comprising both a charged and a neutral sector: 

(3) 

The direction of motion of the fields is fixed by the ex- 
ternal magnetic field. Here we chose for convenience to 
set curvilinear abscissas in such a way that all fields have 
the same chirality (right movers). The periodic neutral 
plasmonic mode for the edge I is 

# (*) = E \f^(% kl e ihx + a[ ki e- ik ne- ha/2 , (4) 

fc;>0 

where a is an ultraviolet cutoff and bosonic creation and 
annihilation operators obey [a ik ,Oj k ,] — SijSkk' with a 
quantized wavevector ki — m^?, m G N. The charged 
zero mode reads 

0?^) = ^rvnix-xi (5) 

with ni the excess number of quasiparticles and xi 
an Hermitian operator conjugate to ri;. The canoni- 
cal commutation relation [xj, ni] — iSji together with 
[^(x), 4>\ {x')\ = ensure that the field on each edge 
satisfies 

[4>i{x),(j)i{x')] = imysgn(x - x'). (6) 

Charge v quasiparticle fields are now defined through ex- 
ponentiation, 

ipi (x) = -jL= e i4 " {x) e i7rux/c > , (7) 
V27ra 

where the extra phase has been added to preserve the 
correct twisted boundary conditions— for the qp field, 
tpl(x + Ci) — Tjji(x)e l2 ^ nil/ . Equation © guarantees that 
the fields -0; create fractional charge excitations. The 
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same commutation relation also insures that quasiparti- 
cles have fractional statistics 

ifo(*MV) = i>}{x')Mxy™ ssnix - x,) - (8) 

For quasiparticles of different edges, fractional statis- 
tics is introduced with suitable commutation relations 
IXjiXl] = inissgn(wj - wi) with w R = -w L = 1 and 
w A = 0. 

With prescriptions @ and ([5]), Eq. ^ becomes 



H? =E l c nf + J2seial s a l<s , 



(9) 



where E c — irvv/Ci is the topological charge excitation 
energy; for the neutral sector, ei = 2itv/ £i is the plas- 
monic excitation energy. 

The total excess charge on an edge is given by 



h = e 



Ci/2 



2n 



= nie 



(10) 



and is conserved in the absence of tunneling. 

Finally, the limit £l,£r — > oo is taken. In the follow- 
ing, for brevity's sake, we relabel the antidot variables as 
C = Ca, e = £a, E c = E^ and n = ua- 



B. Aharonov-Bohm coupling 

The Hamiltonians Eq. ([2]) are expected to describe 
the system anywhere on a v = 1/p plateau since they 
are based on incompressibility. Nevertheless the antidot 
edge, encircling a finite area, is sensitive to the actual po- 
sition in the plateau through a coupling to the Aharonov- 
Bohm (AB) flux^iM We model this effect with an extra 
magnetic field pointing in the opposite direction with re- 
spect to the background magnetic field ^ The AB vector 
potential along the antidot edge reads \A\ = J where 
$ is the AB flux of the additional magnetic field. The 
Aharonov-Bohm coupling is H AB cx j ■ A and describes 
the coupling of the antidot current density with A. By 
a gauge transformation, it is easy to see 2 ^ that this 



amounts simply to a shift in the energies in H A in Eq. §§§ 
according to E c n 2 — > E c (n — $/$ ) 2 , where $o = hc/\e\ 
is the flux quantum. 



C. Tunnel coupling 



Each xLL supports several excitations other than the 
single quasiparticle Eq. 0, given in general by ip™(x) oc 
exp[im<pi(x)], meN. The electron corresponds to m = 
1/v, while the single quasiparticle fields ipj( x ) are ob- 
tained by setting m = 1. One should therefore con- 
sider all possibilities for tunneling, i.e. all terms like 
tl m )^(x)i>T(x), j = L,R. 



Renormalization group (RG) flow equations have been 
set up for the antidot geometry in Ref. 34. Here we 
remind that the renormalized m-quasiparticle tunneling 
amplitude iron scales as a power law 



t 



(m) _ 



(^ren \ ^ 
a J 



('") 



(11) 



under increase in the unit-cell size a — > a ren . Considering 
e.g. 5 = 1/3, the single-qp tunneling amplitude t qp = till, 

diverges when scaling to lower energies, while the elec- 
ts) 

tron amplitude t e = trek scales to zero. The largest a ren 
attainable is the minimum between the thermal length 
oc v/T and the antidot length £. A crossover tempera- 
ture fee Jo ~ E c = nvv/C exists such that for T > To 

one has /t^ m) oc T™* 9 - 1 , while for T < T the flow 
is cutoff by the energy associated to the finite size of the 
antidot ^.cnV< (m) oc E'f 9 - 1 . In the following we will 
assume k^T <C E c and bare tunneling amplitudes such 
that electron tunneling can be neglected. Indeed this ap- 
pears to be the case in most experimental observations 
where single quasiparticle tunneling is clearly observed. 26 
We therefore only retain the dominant term 



j=L,R 



1 E (*jAfa)M0) + h-c> 

3 =L,R 



(12) 

that represents the single-quasiparticle tunneling be- 
tween the infinite edges and the antidot. Here the ve- 
locity v is introduced to have dimcnsionless tunneling 
amplitudes tj. 

A finite source-drain voltage V is applied between the 
left and right edges, producing a backscattered tunneling 
current I(t) of quasiparticles through the antidot 



/(*)= Q L (t)-Q R (t) /2 



(13) 



with Qj(t) = i[Qj,H] = — ie*v[tjip' A (xj,t)ipj(0, t) — h.c] 
and j = L, R. We will consider asymmetric voltage drop 
with aV and (1 — a)V the voltage drops across the left 
and right barrier and < a < 1 (see Fig. [2^,)). 



III. METHOD 

A. Sequential tunneling rates 

For small tunneling as compared with temperature, 
transport can be safely described within the sequential 
tunneling regime^ 6 - Here, the main ingredients are the 
incoherent tunneling rates Tl^h(E). They are obtained 
from the transition probability between the antidot state 
with n qp at time and the state with n' qp at time t, 
at second order in the tunneling Hamiltonian^l 



) P -W j (r) irE 



(14) 
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FIG. 2: a) Double-barrier system, with the voltage drop 
aV [(1 — q)V)] on the left [right] barrier. The forward 
transitions, corresponding to qp transfers contributing pos- 
itively to the current, are indicated, b) Plot of ^fj(x) in 
Eq. (|19|l as a function of fix for g = 1/3,0.8,1, in units 
\tj :\ 2 (u) c /4n 2 )(f3L0c/2it) 1 ~ 9 /T(g) and normalized to its max- 
imum value, c) Scheme of transport regions in the (V, cp) 
plane. Roman numbers indicate the number of charge states 
involved in the transport. Thin lines signal the onset of tran- 
sitions, where energies E± = 0(see Eq. (|24l) ~). 



where 

W,(t) = -(tf(x,t)fi(x t 0)) + (#(*,())#(*, 0)) (15) 

are the thermal correlation functions for the edge (L,R) 
and antidot (A) plasmonic excitations. They can be cast 
in the standard dissipative formal 

Wi{t) = J doj^^- ^(coswt- l)coth^ -isinojt^j 

(16) 

with the antidot and lead spectral densities 

oo 

J A (w) = gve^T 5(u- se)e- u / u ° (17a) 

s=l 

JlMlo) = 9Loe~ u/ ^ (17b) 

where lu c = v/a is & high energy cut-off and f3 = 1/k^T. 
It is important to observe that this expression holds with 
the assumption that the plasmonic fields (f>f are fully re- 
laxed to thermal equilibrium. Mechanisms that could 
guarantee this assumption include for instance thermal- 
izing interactions with external degrees of freedom. In 
the standard xLL theory g = g = v for all edge fields. 
Here we consider the possibility that g, g ^ v to describe 
renormalization effects. 

In the present paper, we do not enter into a microscop- 
ical derivation of exponents and only assume phcnomcno- 
logically that renormalization of exponents takes place. 
The functional form of the quasiparticle correlators is 
thus preserved but is governed by a parameter g = vF . 
The explicit value of F will depend on the details of in- 
teraction, and we will consider it as a parameter. 



Several mechanisms have been proposed to account for 
renormalization. A striped phase is analysed in Ref. [hH 
where interaction is assumed between a v = 1/p xLL 
and phonon-like excitations carrying no net electric cur- 
rent. No tunneling of charge takes place between the %LL 
and the extra phonon modes, nor between the latter and 
the contacts. This assumptions are relevant in the case of 
tunneling between two edge states across a Hall liquid, as 
in the weak backscattering limit of the point-contact ge- 
ometry or of our edge-antidot-edge setup. The only role 
of the additional modes is a modification of the backscat- 
tering dynamics between two chiral edges, leading to a 
renormalized exponent g = i/F, where F > 1 is a func- 
tion of the coupling strength and of the phonon sound 
velocity. 

Another possibility for tunneling renormalization re- 
lies on density-density interactions between chiral edge 
states in a split Hall bar geometry, with parameters g\ 
and gi describing respectively coupling across the con- 
striction and across the Hall bar^ Correlation functions 
for quasiparticles can be calculated exactly in this model 
and give rise to an interedge tunneling current governed 
by g = i>F, where now F depends on g\ and gi. 

Other proposals have been put forward and point to- 
ward more profound modifications of standard chiral Lut- 
tinger Liquids. Interactions with phonons have been con- 
sidered, causing the bosonic field <f> to split into several 
normal modesj 15 i 16 Edge reconstruction offers a further 
scenario, with a non-monotonic density profile of the 
two-dimensional Hall droplet near the edge essentially 
induced by electrostatics! 10 i n ' 19 i 20 ' 21 i 22 

We now return to the expression of the rates Eq. (I14|) . 
It is well known within the bosonized description of edge 
states and reads^^^^ 

T j (E)= £ w slj {E-se), (18) 

S — — QQ 

where 

. (x) _ u ,2 ^ f^A 1 - 9 \T(g/2 + if3x/2n)\ 2 0x/2 
W - \h\ (27r)2 J r{g) e 

(19) 

with r(x) the Euler Gamma function. The sum in 
Eq. (|18[) represents the contribution of plasmons, with 
wcig ht factors w s M At T = they are 

»■ = w(£)''-" Me ">- (20 ' 

with 0(s) the Heaviside step function. 

It is apparent that the rates have contributions from 
both the edge and the antidot correlators. Their func- 
tional behaviour is greatly influenced by the value of the 
lead Luttinger parameter g. In particular, for g < 1 they 
present a non-monotonic behaviour as a function of en- 
ergy (Fig. [lb)). 

Note also that the fractional charge e* is solely deter- 
mined by v and is thus separated from the dynamical 
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behaviour governed by g: it is this separation that allows 
to find independent signatures of v and g. 

We now specify the temperature regime where our se- 
quential tunneling picture holds. 36 A higher limit is set by 
the condition k^T <C E Cl necessary to have a well-defined 
number of qp in the antidot against thermal fluctuations. 
We then observe that the rates scale at low temperatures 
as T 9_1 at energy E = 0. Hence, the linear conductance 
maximum is G max = CGo(/3ui c ) 2 ~ 9 , with Go = e 2 /h and 
G a constant parameter M- For g < 1, G max increases with 
decreasing temperature, implying that for extremely low 
temperatures transport is better described in the oppo- 
site regime of weak electron tunneling^ Sequential tun- 
neling approximation holds if G max <C Go- This implies 
T 3> T m ; n , with T m i n that can be extracted from a knowl- 
edge of typical measurements of G ma x/Go at fixed tem- 
perature. For instance, from experiments^ 6 - performed at 
T ps 10 mK with G max /G ~ 1(T 2 and v = 1/3, assum- 
ing g = g = v we estimate T m ; n pa 1 mK and E c ~ 120 
mK, so that the range between lower and higher limits 
Tmin *C T <C E c /k B spans two full decades in temper- 
ature. The validity range can be extended for systems 
with smaller antidot size (increased E c ) and weaker edge- 
antidot couplings (decreased G). 



B. Moments 

Hereafter, we will introduce higher current moments 
as a tool to determine the x^L exponent and the carrier 
charge. In particular, we will consider the current and 
the p-th normalized current cumulant,— for p — 2, 3 



(I)p 



\eP-i(I)\- 



(21) 



Here, (I) p is the p-th irreducible current moment and 
(7) is the stationary current. They can be expressed- - 
in terms of the irreducible moments of the number 
N of charges e* transmitted in the time r: (I) p = 
lim T _»oo \e*\ p (N) p /r. Observing that (I) = (J)i one has 



p-i 



(N)p 
(N)i 



(22) 



Fano factor and normalized skewness correspond to ^2,3 
respectively. They are expressed as a product of two 
contributions: one coming from the charge of the carrier 
and the other from the statistics of the transport process, 
given in terms of particle number irreducible moments. 

Our task will be to find conditions where (N) p /(N)x 
assumes universal values, independently from e.g. the 
tunneling amplitudes tj and the asymmetry rj = 
\tn\ 2 /\tL\ 2 - We put particular emphasis on the fact that 
such universality must hold independently of the Lut- 
tinger parameter g, since renormalization processes are 
almost invariably present and not controllable. We define 
therefore special the conditions in the parameter space 
where (N) p /(N)i take universal values. Note that the 



statistics of a transport process is identified by all its cu- 
mulants, and therefore all of them should be required to 
be universal to identify special regimes. Here, we will 
adopt only the minimal comparison of the second and 
third moment that are more accessible in experiments. 



C. Master equation 

The detailed analysis of DC current and £2,3 is ob- 
tained from the cumulant generating function calculated 
in the markovian master equation framework^ in the se- 
quential regime. The occupation probability of a fixed 
number n of antidot quasiparticles is 3 -^ 



dt 

n' 3 

(23) 

The energies in the tunneling rates are the differences 
between the antidot and edge j energies before and after 
the tunneling event, 



j^jn — >n+l 



e*V/2 + 2E c ((p-n - 1/2) = El 



E n R +1 ^ n = e*V/2 - 2E c {tp - n - 1/2) = El, (24) 

with if = $/$ + (a— l/2)e*V/2E c . These forward tran- 
sitions n — > n + 1, n + 1 — ► n take place on the left and 
right barrier respectively (see Fig. [2^,)). The correspond- 
ing backward energies obey 



+ 71+1 



Lin — m+1 

— L 

TTin+l — >n 



Y 3 {E)/\t 3 \ 2 , we define T'| 



With T(E) 

the corresponding backward rate T 



—n 

± 



(25) 

= r (£75) and 

r 



De- 



tailed balance states that r ± = e~P E ±T±. 

The conditions E± = grid the (V, ip) plane into dia- 
monds according to the scheme in Fig. [5fc). For the sake 
of clarity, in the following we consider a symmetric volt- 
age drop a = 1/2 when ip depends on the dimensionless 
magnetic flux <!>/<I>o only. However, all results are cast in 
generally valid form. 



IV. RESULTS 

This Section is organized as follows: we describe cur- 
rent, Fano factor and skewness assuming g = g, and we 
discuss various limits in the temperature, voltage and 
flux range, referring to the diamonds in Fig.[2j;). 



A. Current 

A great deal of information on the interaction parame- 
ter can be gathered from the analysis of the DC current, 
although it is not possible to extract signatures of the 
fractional charge. In this subsection we consider g = v, 
despite the generalization g ^ v is straightforward. 
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FIG. 3: Plot of the average current as a function of the 
magnetic flux ip for a) g = v = 1, b) g = v = 1/3 at 
feT = 0.1-Ec- A zoom on the shape of the resonance peak is 
plotted for c) g = v = 1 and d) g — u = 1/3 with different 
temperatures &bT = 0.01, 0.005, 0.001i5 c . Other parameters: 
e*V — 0.05E C , r\ = 1. Curves in panels a) and b) are nor- 
malized to the on-resonance value, while curves in panels c) 
and d) are normalized to the maximum value of the current 
at k B T = 0ME c (solid line). 



1. Magnetic flux dependence 

In the low- voltage regime e*V < 2E C (diamonds I and 
II in Fig.[5J;)) a simple analytical form for the current can 
be found; here, only two adjacent antidot charge states 
are connected through open rates and charge transport 
takes place in a sequence of n — > n + 1 — * n processes. 

In regions I transport is exponentially suppressed be- 
cause of the finite energy E c ; in regions I, II current at 
fixed source-drain voltage oscillates as a function of (p 
with a periodicity of one flux quantum $o for & n V v an d 
g, in accordance with gauge invariance4£ This is rep- 
resented in Fig. [3^) and b) obtained from a numerical 
solution of the master equation. 

Because of periodicity in ip, we start at n = 0. The 
current is 



(I)=v\e*\\t L \- 



,r°r°/_( e *y) 



(26) 



where 



r tot = T° + f + (E° + ) 

f ± (x) = l±e~P x . 



V^_f + (E°_), 



(27) 



If temperature is lowered until k^T <C e*V, then a 
qualitative difference appears in the shape of the reso- 
nance peaks as a function of g. When g < 1, due to 



the non-monotonic rates the current develops two side 
peaks across the resonance and develops a minimum on 
resonance (Fig. [3ji)). On the other hand, for g = 1 the 
rates Eq. (jTSJ) are Fermi functions and the current peak 
displays no structure regardless of the ratio e*V/kBT. 
The only effect is that for lower T plateaus appear with 
increasing width (Fig. Et)). 

It is worth mentioning that the double-peak structure 
for g < 1 was already established in resonant qp tunnel- 
ing through localised impurities,— although in our case 
the parameter that tunes the resonance is the magnetic 
flux. 

If bias voltage is increased, one enters regimes where 
more charge states participate to transport. For voltages 
2F C < e*V < AE C one needs to consider two or three 
charge states. As an example of regions III, we discuss 
the diamond where the rates rj_,r ± , p — n, n + 1, are 
open to transport and the states with n, n + 1 and n + 2 
qp are involved. Again, periodicity allows to choose the 
starting point n = 0. Furthermore, in our temperature 

regime only two backward rates can be retained, T_ and 



T+ (Fig. II left panel). One finds: 

2 i*rlr? + r° + r°_rl 



(I)=r,\e*\\t L \< 



r°(i* +r t 1 ) + ?7 ror 



1 ! 
t 



(28) 



r, . Here also 



where T° t = T° + + r]T_ and rj = rF* 
double-peak structures appear at fixed bias voltage as 
a function <p } although they have a more complicated 
behaviour. Indeed, close to the diamond onset, e*V ~ 
2E C , the current for g < 1 develops a broad minimum 
around if — 1 for large temperature, then a maximum 
in an intermediate range, and finally a minimum again 
with sharp side-peaks at (3E C 1 (Fig. [4j right panel 
b)). No structure appears for g = 1, where decreasing 
the temperature has the only effect of shrinking the peak 
width and causing a larger plateau (Fig. [4] right panel 
a)). 



1.15 
1.10 
1.05 
1.00 





e*V/E c 



FIG. 4: Left panel: scheme of the retained rates in the ap- 
proximation for the three-state region. Right panel: a) cur- 
rent versus ip at g = v — 1 for f)E c = 10, 25, 200, with solid, 
dashed, dotted curves respectively. The curves are at fixed 
e*V/E c = 2.2 and r\ = 1, and are normalized to the value 
of the lowest temperature curve at <p = 1. b) the same with 
g = u= 1/3. 
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2. Bias voltage dependence 

A plot of the current is presented in Fig. [5] as a function 
of the source-drain voltage for symmetric barriers. Again 
the rate behaviour for g = 1 or g < 1 changes qualita- 
tively the current. While for g = 1 the current reflects the 
Fermi-liquid nature of the leads and increases in steps, 
non-monotonic features appear for g = 1/3 and lead to 
negative differential conductance (Fig. [5]b)). Two fea- 
tures are most remarkable in this sense: the on-resonance 
peak (curve at tp = 0.5) for small voltage and the off- 
resonance peak (curve at tp = 0) at e*V « 6E C , indi- 
cated with arrows in Fig. Ob). The first peak develops 
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FIG. 5: Current as a function of voltage for g = v = 
1 (top) and g = v = 1/3 (bottom). The lines corre- 
spond to tp = 0,0.2,0.5 and are shifted of 0.0015/ (0.15J ) 
with respect to the curve at tp = 0. Other parameters 
are r\ = 1 and feT = 0.1J? C . Current units are Jo = 
\e*\Lo c \t L \ 2 {PLo c f~H2^Y-- A /T{g). 



in region II and can be described with Eq. 
tp = 1/2 and r\ = 1, one finds 



r° 



Setting 



(29) 



This result allows for a determination of the Luttinger 
parameter g from the power-law behaviour of the rate, 
r° oc (Pe*V/4n) 9 ~\ valid for /3e*V > 1. 

The second feature finds a natural explanation in terms 
of antidot plasmonic excitations. These in fact enter into 
play for e*V > e = 2E c /v. This new mode increases 
non-monotonically the tunneling rate, leading to peaks 
in current. 

To ascribe the onset of the peak to plasmons, we 
focus on the above example with g = v = 1/3 and 



e = 2E c /v = 6E C . Here, one has to consider five charge 
states, n, . . . , n + 4 and the rates T± , p = n, . . . , n + 3. In 
this regime it is possible to obtain an approximate form 
of the current neglecting all backward rates except those 

Q 3 

opening at the border of region V (r_,T + , see Fig. El 
inset). We work at tp = 2, choosing n = as a reference. 
Observing that the large-energy behaviour of the forward 
rates is oc (f]E) a ~ 1 , we can process the rates according 
to the magnitude of their argument in region V: further- 
more, for n = l,T l + = rl"' with I = 0,1, 2, 3. Thus we 
set Tj_ = Tj_ ~ /, with / a temperature dependent con- 
stant, since these rates are in their slow-decaying regime. 
On the opposite, the full energy dependence of = r°_ 
is taken into account. In terms of the 'plus' rates only, 
the current is given by 



f r 

(I) = 2\e*\ J - 



r° + (r 3 + + /) 



3r° + 2r3 



3Ft 



(30) 



Finally, we set T5_ 



/ + 5T p , where ST p is the first 

(I) 
0.5 
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3 




/<+ / n; 


2 








X 


1 






2 




v e*V J 




2 4 







5.25 5.5 



5.75 
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FIG. 6: Current in the IV-V state regime versus bias voltage 
at g — v — 1/3. Solid line, numerical result; dot-dashed line, 
approximate solution in Eq. (|30[) with 8F P = gT+; dashed 
line, approximate solution without plasmons (ST P = 0). The 
curves are at tp — 2 (equivalent to tp — in Fig. [5] be- 
cause of the periodicity in magnetic flux), feT = 0.175c, 
rj — 1 and units are as in Fig. [5] The best fit parameter 
is / = 0.19Jo/|e*|. Inset: scheme of retained rates in the 
approximation for region V. 

plasmon contribution to the rate that shows up as a sec- 
ondary peak at energy E = e. We now compare the 
current with no plasmons (ST P = 0) with ST P = gT^_ 
obtained from Eqs. ([T5]) and (f2T)]) with relative weight 
ivi/wo = r(<? + 1)/T(g) — g . Figure [6] shows a com- 
parison between the current with these approximations 
and the result obtained through a numerical solution of 
the master equation. It is apparent that neglecting the 
plasmon completely misses the peak onset that is instead 
well-captured by out latter approximation. This testifies 
that plasmonic excitations play indeed the crucial role in 
current enhancement. 
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B. Current moments 

Current moments are now discussed. 25 According to 
Eq. (22), these two quantities are the natural observables 
to look at in order to measure the quasiparticle charge 
e*. In the antidot geometry, different conditions can be 
found where this measurement is special in the sense dis- 
cussed above. We remind that for a poissonian process, 
as with weak backscattering current in a point contact, 
k 2 =e*/e = u, k 3 = (e*/e) 2 . 



1. Few-state regime: e*V < 2E C 

For the sequential tunneling master equation Eq. ([23]) 
restricted to the two-state regime (two charge states) , an 
exact analytical treatment is possible. 

A known formula 38,44 for the Fano factor is obtained 



coth 



(3e*V 



2 V 



r° + r°_f_(e*v) 



r 2 

1 tot 



(31) 



while for the skewness we find 

2 



fa = - 



1 — 6i] 



T° + T°_f + (e*V) 



r 2 

1 tot 



,rfr° 2 / 2 (eT) 

-12ry 2 + 



r 4 

1 tot 



(32) 



with the functions r to t, f± defined in Eq. ([27)1 . The pre- 
vious equations are an example of the intertwining be- 
tween charge and process statistics that does not allow, 
in general, for unequivocal conclusions on e* . We analyse 
now the behaviour varying the ratio e*V " jk^T looking for 
special conditions. 

Thermal limit: e*V -C k^T . The Fano factor is indepen- 
dent of the charge fractionalization, 



, fc B T 
' eV ' 



(33) 



reflecting the fluctuation-dissipation theorem. On the 
contrary, the normalized skewness that measures the fluc- 
tuation asymmetry induced by the current depends on 
the carrier charge e* = ve. Indeed, for low voltages 
V — > + one has 



fa 



V 



1 



(l + r/) 2 Cosh 2 (/?£ c (<p-l/2)) 



(34) 

Note that the tp dependence can be used to extract r\ and 
e*/e independently from g. 

Shot limit: k^T <C e*V. In the blockade regions I with 
|/3-E±| ^> 1, one has fa = v and fa = v 2 . In this case 
the statistics of the transport process is poissonian: the 
transport through the antidot is almost completely sup- 
pressed, (I) w 0, and the residual current is generated 
only by a thermally activated tunneling that is com- 
pletely uncorrelated. So regions I constitute an example 
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FIG. 7: Fano factor (left) and skewness (right) as a function 
of (p. Panels a) and b) r\ = 1, panels c) and d) r\ = 3. Param- 
eters: e*V = 0.LE c , k B T = 0.004_E C and 3 = 1/5 (solid line), 
g — 1/3 (dot-dashed line), g — 1/2 (dotted line). Horizontal 
thin lines indicate the universal limits at 1/2 at 1/4. 



of special regime. We consider now the two-state regime 
(II) for (3E± 3> 1. For fractional edges g < 1, fa^ have 
a particular functional dependence on ip. We find that 
they both develop a minimum and that for not too strong 
asymmetries the absolute values of the minima are 



2 n- 



(35) 



These minimal values do not depend neither on g, as the 
comparison of solid (g — 1/5), dotted (g — 1/3) and 
dashed (g = 1/2) curves in Fig. [7] confirms, nor on r\. For 
Fermi liquid edges g = 1, we have 



fa 
fa 



v 2 [l - 6*7(1 - 



V)/(i 



■v) A 



(36) 



independently from tp. Here, fa and fa assume their 
minimal values v /2 and ^ 2 /4 in the symmetric case r) = 1. 
In this conditions we have the strongest anticorrelation 
that is signalled by a marked sub-poissonian statistics. 

We can conclude that in the two-state regime, in the 
shot limit, the values of the minima for fa^ obtained 
varying rj, ip correspond to a special condition where the 
system shows the same universal sub-poissonian statis- 
tics for any g < 1. This represents a means of testing 
fractional charge outside poissonian conditions and in- 
sensitive to renormalizations of the Luttinger parameter. 

In the intermediate regime e*V w k^T, fa^ depend 
more strongly on the parameter g, and the interplay 
of these two energy scales prevents the onset of special 
regimes. 
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FIG. 8: Panel a): Fano factor k^/v at u = g = 1/3, ksT = 
0.02_E C , vs. source-drain voltage and ip. Panel b): normalized 
skewness k^/v 2 with the same parameters. Both moments 
are plotted for symmetric barriers rj = 1 (top) and strong 
asymmetry n = 10 (bottom). 



J2. Many-state regime: e*V > 2_B C 

We study now higher voltages e*V > 2E C where the 
renormalized interaction parameter g has a prominent 
role. For this purpose we consider the behaviour of both 
the Fano factor and the normalized skewness. In Fig. [5] 
a density plot of ki and ks for v = j = 1/3 as a func- 
tion of magnetic flux and source-drain voltage is shown 
for different asymmetries. First of all, in region I we 
recover poissonian statistics: k-2 = v, k% = v 2 (middle 
grey) for any rj. Outside this region, light grey zones 
represent super-poissonian Fano factor and normalized 
skewness (&2 > v and k^ > v 2 ), while dark grey re- 
gions represent sub-poissonian behaviour (&2 < v and 
^3 < v 2 )• Figure [8] shows that super-poissonian regions 
are possible with g < 1. For g > 1, we always have 
sub-poissonian behaviour in accordance with previous re- 
sults (not shown).— In presence of asymmetry, super- 
poissonian regions increase. We note that Fano and skew- 
ness present concurrent super/sub-poissonian behaviour, 
and that the maximal values of the skewness in the super- 
poissonian regions are stronger. Given this similarity, in 
the following we will discuss the Fano factor only. 
Three-state region. In region III, a tractable analytical 
formula for the Fano factor can be derived under the 
same assumptions made for the current Eq. (|28|) : one 
has k2jv = 1 — 1r\ 8k2, with 



Sk 2 = 



rfr^ri 



T°r° 



-r°_rV (r?-rj) forl-r°; 



[vT°-ri 



(37) 



with r£,Tj in Eq. {28]). We note that in order to have 
super-poissonian noise a fractional g < 1 is necessary, 
with additional conditions on the asymmetry. Indeed, 
setting r) = 1 in Eq. (J37J) in the limit (3E\,f3E Q _ > 1 
yields 8k2 > for any g. On the other side, setting g = 1 
gives <5&2 = 2rj / (rj 2 + rj + 1) 2 > 0. It appears that positive 
correlations are induced by an interplay of rj and g. 
Five-state region. Finally, interesting effects take place 
in the five-state regime (V) for v — 1/3. Here, a strongly 
super-poissonian Fano factor appears along the diamond 
lines for rj = 1 and disappears for large asymmetries. 
An investigation of this effect can be performed with 
the same methods used to obtain Fig. [6] The rate en- 
hancement due to the onset of the plasmonic collective 
excitations can again be shown to be responsible of the 
super-poissonian behaviour at small asymmetries. 



V. CONCLUSIONS 

In conclusion, we have analysed transport of quasipar- 
ticles through an antidot coupled with edge states in the 
fractional Hall regime. The model of a finite size chiral 
Luttinger liquid with periodic boundary conditions has 
been reviewed and cast in a suitable form for calculations 
of higher current moments through a master equation ap- 
proach in the sequential regime. We have also allowed for 
the possibility of a phcnomcnological renormalization of 
the interaction parameter. 

We have found that independent information on in- 
teraction renormalization and fractional charge can be 
extracted from tunneling current and its moments, noise 
and skewness. For current, remarkable qualitative dif- 
ferences as a function of the Luttinger parameter appear 
in the shot limit of resonance peaks and in the three- 
state regime. A quantitative determination of g is fur- 
thermore possible through the power-law behaviour of 
on-resonance current versus voltage at low temperatures. 
Current moments also depend strongly on the interaction 
parameter: in particular, super-poissonian behaviour is 
never found for g > 1 regardless of asymmetry. On the 
other hand, we have identified special regimes in the one- 
and two-state regions where a comparison of Fano factor 
and normalized skewness realizes an unambiguous charge 
determination procedure, insensitive to renormalizations 
of g. Finally, signatures of plasmonic excitations are in- 
dicated in the large bias voltage regime. 

Confirmation of such novel results appears to be within 
reach: plugging estimated parameters from present 
experiments^ into our results gives currents in the range 
of 0.5 5 pA for .9 = 1/3 and T = 10 mK. At the same 
time, recent accomplishments in measurement techniques 
applied to electron counting open the possibility for fea- 
sible noise and skewness determination, even in systems 
with very low current and noise levels 4^ 

Furthermore, litographic approaches have been devel- 
oped that allow for antidot radii sensibly smaller than 
the 300 nm in the original experiments^ Thus, the en- 
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ergy scale E c associated with the antidot finite size can 
be raised to the order of some hundreds of mK, and more 
easily attainable temperature regimes can be compatible 
with the requirements of our approximations. 

We believe that our results help clarify some issues re- 
lated to transport in the Hall regime, namely features 
like fractional charges and nonpoissonian correlations. 
Such a thorough analysis of a Hall antidot device is a 
necessary building block, especially in view of possible 



applications of similar systems to topological quantum 
computation^ 
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